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1 Introduction 

Classical theory of wave turbulence was developed for description of statis- 
tic properties of weakly nonlinear waves in infinite domains. Importance of 
resonant interactions among these waves was first pointed out in 1961 by 

0. Phillips pQ and this led to construction of the first kinetic equation |2j. 
Later on it was established by many researchers that effects of finite length 
are not described by the classical theory and kinetic equations, and 20 years 
later, in 1981, O.Phillips wrote that "new physics, new mathematics and new 
intuition is required" |3j in order to understand energetic behavior of these 
systems. Discrete effects appear when waves "notice" boundary conditions, 

1. e. wave lengths are comparable to the sizes of interaction domain, and 
therefore wave numbers are integers. First theoretical result on the behavior 
of discrete waves systems was published in 1990 in [3J where main distinction 
between continuous and discrete wave systems was established: stochastic 
interactions between all waves (infinite domain) versus interactions within 
small independent groups of waves and existence of non-interacting waves 
(finite domain). Theory of discrete wave systems was then presented in jS], 

[Zj thus giving a rise to the following qualitative picture of wave turbu- 
lence: short waves are described by power energy spectra and kinetic equa- 
tions, long waves are described by Clipping method and dynamic equations. 
Numerous papers of last few years showed that this qualitative picture has to 
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be modified because some discrete effects (frozen turbulence, mesoscopic tur- 
bulence, etc.) are also observable in the short-waves' part of the wave spectra 
(|H], [H] ^Dj and others). Discovery of intermittent patterns in that part of 
wave spectra which is supposed to be described by a kinetic equation only, 
gave us the reason to revisit the very foundations of wave turbulence theory. 
It turned out that its basic theoretical ground - KAM theorem - is proven 
not for all short waves, a countable number of short waves with rationally 
connected frequencies is excluded from consideration leaving some "gaps" 
in the wave spectrum. These gaps correspond to the waves whose ratio of 
frequencies is algebraic number of degree < 2, for instance, a rational num- 
ber which is usually understood as necessity for the corresponding dispersion 
function to be a rational function. A model of laminated turbulence was re- 
cently presented ^T] which consists of two layers, discrete one (in the whole 
range of wave numbers) and continuous one (in the short-waves' part of the 
wave spectrum). The main fact allowing to construct this model is following: 
frequencies of discrete resonantly interacting waves are rationally connected 
though the dispersion function itself can be highly irrational. Thus, 
the model of laminated wave turbulence "fills" the gaps left by KAM theory 
and explains in particular coexistence of power energy spectra and coherent 
structures in the short-waves' part of wave spectrum. 

Model of laminated turbulence brings our attention to the new question: 
how to compute discrete resonances in large computational domain? Indeed, 
say for 4-waves interactions of 2-dimensional gravity waves, the resonance 
conditions can be regarded in the form 

y/ki + \fk~i = yfa + y/k±, ki + k 2 = k 3 + k A 

where ki = (m^nj), \/i = 1,2,3,4 and fej = \ki\ = y/mf + It means that 
in a finite but big enough domain of wave numbers, say \m\, \n\ < 1000, direct 
approach leads to necessity to perform computations with integers of order 
10 6 . These computations in a substantially smaller domain \m\, \n\ < 128 
took 3 days with Pentium-4 (|12j). 

The main goal of this paper is to show that constructions of indepen- 
dent classes of resonantly interacting waves jH] can be re-formulated as a 
pure mathematical procedure which allows to reduce drastically computa- 
tion time for this sort of equations. Gravity waves are taken as our main 
example. Some necessary results on discrete resonances of gravity waves are 
given in jH] but formulations there are too formal and partly inaccurate. We 
re-formulate all necessary results on construction of classes for these waves in 
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Sec. 2 and give some illustrative examples. In Sec. 3 we describe the compu- 
tation algorithm and some preliminary results of numerical simulation. Brief 
discussion is given in Sec. 4. 

2 Construction of classes 

In this section we regard 4-waves resonances of 2-dimensional gravity waves 

in the form 

\fk~i + \ffa = ^ffa + v^4 (1) 

fa + fa = fa + fa. (2) 

To construct independent classes of resonantly interacting we need following 
definition. 

Definition Let k = (m, n) be a vector with integer coordinates, m,n G Z. 
Represent the square root of the norm \k\ as 

(m 2 + n 2 ) 1/4 = 7g 1/4 , 7,9 G N 

and q does not contain fourth degrees: 

q = pT---pT, ate {1,2,3}, 

with all different prime p^. Number q is called index of a vector k = (to, n) 
and number 7 is called its weight. Class Cl q of index q is set of all vectors 
with this index. 

For instance, vectors (1,3) and (12,4) belong to the same class Clio of 
index 10, with weights 7 = 1 and 7 = 2 correspondingly, while vector (2, 1) 
lies in Cl$ of index 5, with weight 7 = 1. 

Lemma Let vectors fa, fa, ^3, k^ construct a solution of (|T|). then only two 
cases are possible: 

Case 1: all vectors belong to the same class, 

Case 2: all vectors belong to two different classes Cl qi , Cl q2 in such a way 
that there exist q\ , q% : 

fa, fa E Cl qi and fa, fa € Cl q2 

or 

fa, fa G Cl qi and fa, fa G Cl q2 . 
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In Case 2 all solutions are symmetric, i.e. vectors belonging to each class 
must have the same norm. All asymmetric solutions, if any, are described by 
Case 1. 

The statement of Lemma follows immediately from elementary proper- 
ties of integers and we are not going to the detailed proof here. The general 
idea of the proof is very simple indeed: two different irrational numbers can 
not satisfy any equation with rational coefficients. For instance, equation 
ay/3 + by/5 = c has no solutions for arbitrary rational a, b, c. It means that 
irrationalities corresponding to classes indices in (JJJ have to be rid off in 
order to construct its integer solutions. In general, one to four different ir- 
rationalities can appear in (JTJ), simple consideration show that all cases but 
two described below give no integer solutions. 

In Case 1, all four irrationalities are equal and can be canceled, so that 
(JTJ) is reduced to 

7i<7 1/4 + 72<7 1/4 = 73<? 1/4 + 74<? 1/4 => 7i + 72 = 73 + 74- 

Let us regard as example a couple of asymmetric solutions given in |lUj : 

(-4, 0), (49, 0), (9, 0), (36, 0) and (-20, 15), (-20, -15), (-49, 0), (9, 0). 

An easy check shows that all 8 vectors belong to the same class Cl\ of index 
1, with weights 71 = 2,72 = 7, 73 = 3,74 = 6 and 71 = 5,72 = 5, 73 = 
7, 74 = 3 correspondingly, so that condition 71 + 72 = 73 + 74 is full-filled in 
both cases. It is important to understand that 4 vectors constructing some 
asymmetric solution must belong to the same class but not necessarily to Cl\ 
as in the solutions above. In the next section some asymmetric solutions are 
presented which lie in other classes. 

In Case 2, the irrationalities with corresponding coefficients in front of 
them must be pair-wise equal to be canceled, for instance 

1/4 1/4 1/4 1/4 , 
7i0i + 7202 = 739i + 74<? 2 7i = 73 and 72 = 74. 

Notice that Lemma gives only necessary condition of the existence of 
a solution. This means that to find a solution we have first to construct 
classes, then find solutions within the corresponding classes and check linear 
conditions afterwards. Obviously, (JU,© have infinitely many solutions - for 
instance, quartets 

ki = (a, b), k 2 = (c, d), k 3 = (a, b), k 4 = (c, d) (3) 
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and 



k\ = (a, b), k 2 = (c, a — b + c) , fc 3 = (b, a), k^ = (a — b + c, c) (4) 

with arbitrary integers a, b, c, d give its solution, as well as any proportional 
quartets corresponding to multiplication of all wave-numbers on the same 
integer. Less trivial example of "tridents" 

fci = (a, 0), k 2 = (-6, 0), k 3 = (c, d), k 4 = (c, -d) (5) 

possess two-parametric series of solutions [TQJ : 

a = (s 2 + t 2 + st) 2 , 6 = (s 2 + t 2 - st) 2 , c = 2st(s 2 + t 2 ), d = s 4 - t 4 (6) 

with arbitrary integer s, t. This series gives solutions of both (HJ,© though 
perhaps not all of them. Parametrization (JBJ) is constructed in such a way 
that norms of all four vectors are full squares, i.e. again vectors hi, k 2 , ks, k± 
belong to the same class Cl\ of index 1, with weights 

7i = s 2 + t 2 + st, 72 = s 2 + t 2 - st, 7 3 = 74 = s 2 + t 2 , 

providing 7i + 72 = 73 + 74 for any parameters s,t. 

There is no known way to construct general solution of (0),(j2J) and even 
a construction of some particular solutions' series demands a lot of skillful 
work and a bit of luck. On the other hand, this hard work is mostly not 
well-paid because coming back to physical problem setting, we usually need 
all solutions in some spectral domain, not only those which could be nicely 
parameterized. In the next section we use Lemma to construct fast computer 
algorithm for finding solutions of (fljl - P]) which is a challenging numerical 
problem. 

3 Scheme of numerical algorithm 

As it was shown in previous section, equations (JU,© have infinitely many 
solutions. But not all of them are interesting from physical point of view. 
We follow [E] in classification (written out up to the cyclic change of indices) 
of all possible solutions according to their role in energy transfer : (I) trivial 
resonances ki = k 3 , k 2 = A; 4 which do not redistribute the energy, an exam- 
ple is given by ©; (II) symmetrical resonances = \ks\, \k 2 \ = \k^\ which 
do not generate new wavelengths and therefore do not transfer the energy 
through the scales, an example is given by (J2J; (III) asymmetrical resonances 



5 



playing major role in energy transfer, an example is given by ©. In this pa- 
per we are primary interested in asymmetric solutions which means in terms 
of Lemma that all vectors belong to the same class (Case 1). Not going into 
all programming details we present here just underlying main ideas of our 
algorithm and first results of computer simulations. 

Step 1. Create first auxiliary array of primes A p in the given spectral do- 
main p < D using Eratosthenes' Sieve procedure which eliminates composite 
numbers from the list of natural numbers and can be briefly described as 
follows: 

• Create an array Ad containing all integers l,...,d, write number "1" 
into array A p , mark "1" in the array Ad as used element; 

• Introduce cycle variable i — 1, Vd; 

• Find first number in Ad which is greater than 1 and is not marked yet 
as a composite, denote it as r and mark all numbers 2r, 3r, 4r, ... as 
composite (at the first step r = 2, at the second r = 3 and so on); 

• write number r into array of primes A v \ 

• put i = r and repeat the procedure. 

Step 2. Create second auxiliary array A q of possible indexes constructed 
of the primes obtained at the Step 1. Indexes are computed directly by 
formula 

Q = pT---pT, a, €{1,2,3}, 

in the domain q < 2 l ^\fd < 1.19V"- Say, for d = 1000, we have q < 37.63, 
i.e. it is enough to compute for q < 37. 

Step 3. Find all integer solutions of the linear equation 

7i + 72 = 73 + 74 
in the domain ^ < 2 l l i \fd 1 Vz = 1, ..4. 

Remark. Further it will be necessary to check when expressions under 
the radicals 7 4 g can be represented as a sum of two squares. Here one has 
to make use of the well-known Euler theorem: representation of a positive 
number as a sum of two squares which is possible iff each of its prime factors 
of the form At + 3 occurs as an even power. For instance, pair 7 = 3, q = 1 
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has 3 in even degree, expression under the integral is 3 4 • 1 = 81, represen- 
tation as a sum of two squares exists and number 81 = 9 2 + 2 should be 
written into the array A 1S . On the other hand, for pair 7 = 2, q = 3, the 
expression under the radical has factor 3 in the odd power, i.e. number 48 
can not be represented as a sum of two squares and it is not an element of 
the array A 1>q . Notice that all prime factors of 7 have even powers which 
means that only prime factors of q has to be investigated. Consequently, an 
arbitrary presentation of q as a sum of two squares, multiplied by 7 4 , gives 
corresponding presentation of the expression under the radical. 

Step 5. Check Euler theorem for all elements of A q and construct ar- 
ray A q with "allowed" elements only; find all presentations as a sum of 
two squares for each "allowed" element. Obviously, one number can be 
decomposed into sum of two squares in more then one way, for instance 
1105 = 4 2 + 33 2 = 9 2 + 32 2 = 12 2 + 31 2 = 23 2 + 24 2 . 

Remark. Due to Lemma, the search of 2-square-representations for the 
numbers < 2<i 2 is reduced to the numbers < \.l§\fd. In particular, for 
d = 1000 we have to compute these presentations for the numbers in the 
domain < 37 which can be done by direct enumeration. The number S of 
different 2-square-presentations of an integer is proportional to the difference 
between its prime factors of the form At + 1 and At + 3 (Lagrange theorem). 

Step 6. For each q G A q and each its possible 2-square-presentation, 
q = mj+nj,j = 1, 2, S, compute all under-integral-expressions 7 4 (m 2 +n 2 ) 
and check linear conditions on rrij, ny. 

m 1 +m 2 = m 3 + m 4 , 

ni + n 2 = n 3 + n 4 , 

with rrii = li m j, n i — li^ji m 2 = 7l 7 %> n i = It,™ 31 m 3 = it™ 31 n 3 = 
ilhj, m A = 74 m.,-, n 4 = 7!%. 

Summarizing, Lemma allows to reduce domain of integers under consid- 
eration from 2c? 2 to 2 1 ^ 4 \^d while classical procedure (Eratosthenes' Sieve) 
and some known number-theoretical results (Euler and Lagrange theorems) 
allow to further reduce the computation time which is about 4 min. at 
Pentium-4 for spectral domain m^rii < 1000. In this domain we have found 
five asymmetric solutions which are not tridents: 

fci = (495, 90), %2 = (64, 128), k 3 = (359, 118), h = (200, 100), 
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fci = (675,225), k 2 = (64, 192), k 3 = (479,237), k A = (260, 180), 

fci = (810,45), A-2 = (128, 192), k 3 = (598, 117), k 4 = (340, 120), 

fci = (855,360), k 2 = (64,256), k 3 = (599,356), k 4 = (320,260), 

h = (990,180), k 2 = (128,256), k 3 = (718,236), h = (400,200). 

These solutions belong to classes C7 5 , C7 10 , Cl 13 , Cl 17: Cl 20 correspond- 
ingly, with weights 71 = 15, 72 = 8, 73 = 13, 74 = 10 for all solutions. 



4 Brief discussion 

We described here first version of our algorithm in order to show how to use 
Lemma for reducing drastically computation time. At present, this algorithm 
is written only for positive integer wave numbers m^rii G N. Using some 
simple symmetrical considerations, one can adapt it for the case of arbitrary 
mi,rii G Z. Another simple modification can be made in order to find solutions 
for the case 



or for the case of 5-wave interactions and more (of course, Lemma should be 
re-formulated) . 

Moreover, basing on this algorithm, one can develop generic algorithm 
for fast computation of discrete resonances among different types of waves 
with dispersion function being an arbitrary polynomial function of the 
norm of wave vector, = 0(k). Different dispersion functions in case of two 
dimensional waves will have different formula for index (Step 2 of present 

— * 

algorithm); in case of three dimensional waves, i.e. k — \k\ = ym 2 + n 2 + I 2 , 
Lagrange and Euler theorems should be replaced by some known number- 
theoretical results on the decomposition of an integer into the sum of three 
squares. All other algorithmic steps will be the same as above. Development 
of such a generic algorithm for a big class of dispersion functions is our cur- 
rent object of interest. 

Author is very grateful to Sergey Nazarenko for valuable discussions and 
remarks. 
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